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Abstract 

Phase separation is a fairly common physical phenomenon with examples 
including the formation of water droplets from humid air (fog, rain), the 
separation of a crystalline structure from an isotropic material such as a liquid 
or even the formation of the sizzling gas bubbles when a soda can is opened. 

It was recognized long ago (at least on a phenomenological level) that sys- 
tems exhibiting several phases in equilibrium can be described with an appro- 
priate variational principle: the phases arrange themselves in such a way that 
the energy associated with the phase boundaries is minimal. Typically this 
leads to an almost deterministic behavior and the phase boundaries are fairly 
regular. However, when looked at from a microscopic point of view, the sys- 
tem consists of a bunch of erratically moving molecules with relatively strong 
short-range interaction and the simplicity of the above macroscopic descrip- 
tion looks more than miraculous. Indeed, when starting from the molecular 
level, there are many more questions to be asked and understood: which are 
the phases which we will see? why do only those occur? why are the phase 
boundaries sharp? how should we find (define) the energy associated with the 
interfaces? Only then can we ask the question: why does the system minimize 
this energy? 

It is only in the last decade that a mathematically satisfactory understand- 
ing of this phenomenon has been achieved. The main goal of the talk is to 
present the current state of affairs focusing thereby on results obtained in 
joint works with Raphael Cerf. The connection to fields of mathematics other 
than probability theory or statistical mechanics will be highlighted; namely, 
to geometric measure theory and to the calculus of variations. 
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1. Introduction 

Although phase separation is a fairly common physical phenomenon (exam- 
ples will be given further below) its mathematically satisfactory understanding, even 
in the simplest models, has not been achieved until the last decade. In order to 
uncover the mechanism leading to the separation of various phases in an initially 
homogeneous material, in particular to explain why and what kind of phases will 
occur and which shapes they take, one has to work with a microscopic description 
of the system, often at a molecular scale. Materials on such scales however tend 
to behave 'chaotic' and this strongly motivates (if not compels) the use of a prob- 
abilistic approach. In this approach, the system is modeled by randomly moving 
(in equilibrium theory randomly located) particles which interact with each other 
according to some simple, typically short range, mechanism. The goal is then to 
derive the large scale behavior of the system based only on the specification of the 
local interaction. 

The difficulty of the analysis stems from the fact that the interaction, although 
local, might be strong enough (depending on some parameter such as the temper- 
ature) to cause a subtile spatial propagation of stochastic dependence across the 
entire system. As a consequence, one has to leave the familiar realm of classical 
probability theory whose focus has been laid on the large scale effects of randomness 
arising from independent (or weakly dependent) sources. Instead, we have to deal 
with a strongly dependent system and it is exactly this strong dependence which 
causes a highly interesting cooperative behavior of the particles on the macroscopic 
level which can, in certain cases, be observed as phase separation. 

The problem of phase separation and related issues have been a driving force 
behind developing, and a benchmark for testing various new techniques. Postponing 
historical remarks until section 2.1, let me highlight here only those ones which play 
an essential role in our approach [10, 11] achieved in collaboration with Raphael 
Cerf. The basic framework is (abstract) large deviation theory, see e.g. [37], whose 
power contributed substantially (admittedly rather to my own surprize) to the suc- 
cess of this approach. To have sufficient control of the underlying model, in our case 
the Ising-Potts model, we employ spatial renormalization techniques, as developed 
in [34], in conjunction with the Fortuin-Kasteleyn percolation representation of the 
Potts model, [21]. Finally, tools from geometric measure theory a la Cacciopoli 
and De Giorgi will be employed to handle the geometric difficulties associated with 
three or higher dimensions, in a similar fashion as was done in [9]. 

The goal of this article is to present some recent developments in the equilib- 
rium theory of coexisting phases in the framework of the Ising-Potts model. The 
presentation will be based mainly on the results contained in [34, 10, 11]. It will 
include a description of the underlying physical phenomenon, of the corresponding 
mathematical model and its motivation and, following the statement of the main 
results, comments on the proofs will be included. In order to address an audience 
broader than usual, I will try to use as little formalism as possible. 



1.1. Phase separation: examples and phenomenology 
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Perhaps the most common and well known example of phase separation is 
the development of fog and later rain from humid air. When warm humid air 
is cooled down so much that its relative humidity at the new cooler temperature 
would exceed 100% (i.e. it would become over-saturated) the excess amount of 
water precipitates first in the form of very small droplets which we might observe as 
fog. The system at this time is not in equilibrium, rather in a so called metastable 
state. After waiting very long time or simply dropping the temperature further 
down, the droplets grow bigger and ultimately fall to the ground due to gravitation 
in the form of common rain droplets. In this example phase separation occurred 
since from a single homogeneous phase (warm hiimid air) two new phases have 
been formed: a cooler mixture of water and air (note; with 100% relative humidity 
(saturated) at the new lower temperature) plus a certain amount of water, more 
precisely a saturated solution of air in water, in the form of macroscopic droplets. 
In fact, in the absence of gravitation and after very long time, only a huge droplet 
of fluid would levitate in a gas (both saturated solutions of air in water and water 
in air, respectively). The opposite situation (water majority, air minority) may also 
occur. Consider the following familiar example; think of a bottle of champagne 
when opened. Here the change of temperature is replaced by a change of pressure 
but the phenomenon is similar with the roles exchanged; at the new lower pressure 
the liquid is not able to dissolve the same amonunt of carbon dioxide, hence this 
latter precipitates in the form of small bubbles (droplets), etc. 

The phenomenological theory explains this type of phenomenon as follows. At 
any temperature there are saturation densities of air/water and water/air mixtures 
and only saturated solutions will coexist in equilibrium. They also determine the 
volumes of the two coexisting phases. Moreover the phases arrange themselves in a 
way so as to minimize the so-called surface energy, associated with the interface be- 
tween the phases. In fact, it is supposed to exist a (in general direction-dependent) 
scalar quantity r, called the surface tension, whose surface integral along the inter- 
face gives the surface energy. The surface tension, as well as the saturation densities, 
have to be measured experimentally. It is implicitly assumed that the interface is 
'surfacelike' and regular enough so that the integral along the surface makes sense. 
The prediction which can be made is that the shape of the phases in equilibrium is 
just a solution of the variational principle. By the classical isoperimetric theorem, in 
the isotropic case the solution is just a sphere, hence the occurrence of bubbles and 
droplets. In the non isotropic case the corresponding variational problem is called 
the Wulff problem. The solution is known to be explicitly given [38, 13, 19, 20, 36] 
by rescaling appropriately the so called Wulff crystal: 

Wt = {x €W'' ; X- u < r(z/) for all unit vectors u}. 

It is worth noting that the same arguments are used to describe macroscopic crystal 
shapes as well. 

1.2. The mathematical model and the goals of the analysis 

The next step is to find a model which is simple enough to be analyzed by 
rigorous methods yet rich enough to exhibit the phenomenon we want to study. 
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In order to accommodate a multitude of phases we consider a finite number (g) 
possible types of particles (called colors or spins). Physics suggests to choose a 
short range interaction of "ferromagnetic" character which means that particles of 
identical type prefer to stay together and/or they repel particles of different types. 
For simplicity we assume that the interaction distinguishes only between identical 
and different types, otherwise it is invariant under permutation of colors. There is 
a standard model of statistical mechanics, called the (ferromagnetic) q-states Potts 
model, which corresponds exactly to these specifications. We consider the closed unit 
cube f2 e M*^, > 3 (modeling the container of the mixture of particles) overlapped 
by the rescaled integer lattice Zf^ = Z'^/n. We define f7„ = 17 n ZjJ, and denote 
by d^"^n the internal vertex boundary of r2„. At each lattice point x there is a 
unique particle of one of the types 1, 2, q. The energy H{a) of a configuration 
<^ = i<^x)xen„ can be chosen to be the mmrber of nearest neighbor pairs of different 
types of particles corresponding to nearest neighbor repulsion. According to the 
Gibbs formula, the probability of observing a configuration a is proportional to 
g-/3ff(o-)^ where /? = 1/T is the 'inverse temperature' which adjusts the interaction 
strength. (High T = large disorder = relatively small interaction, etc.) Note that 
we use a static description of the equilibrium system. It corresponds to a snap-shot 
of the system at a given time and the task is to understand the 'typical' picture we 
will see. 

A restricted ensemble is a collection of certain feasible configurations. For 

instance, in the situation of the water/air mixture every configuration with a fixed 
number of water and air particles is possible, and this collection forms our restricted 
ensemble. It turns out that the direct study of this particular ensemble is extremely 
difhciilt and it is a crucial idea (discovered long ago) to go over to a larger, more 
natural ensemble, namely to that without any restrictions on the particle numbers. 
Then, the restricted system can be regarded as a very rare event =: G„ in the 
large ensemble and conditional probabilities can be used to describe the restricted 
system. The events Gn are often in the large deviations regime, and it is from here 
that large deviations theory enters the analysis in an essential way. 

The unrestricted system is usually referred to as the Potts model with free 
boundary conditions. In the case q = 2, it is equivalent to the classical Ising model. 
The Gibbs formula and the energy uniquely determines the probability measure 
in this (and in every) ensemble. We can introduce mixed boundary conditions as 
follows. Divide the boundary F = ^^l of the 'container' into q + 1 parts indexed by 
F*^, F^, ....F'. The parts can be fairly general but the {d — l)-dimensional Hausdorff 
measure of their relative boundaries has to be zero. We set for n e N and i = 
0,...,q, 

= {xGr„; doo(x,F*)<l/n and Vi<i, doo(a:,FJ)>l/n} i = 0,...,q 

where doo denotes the distance corresponding to the max-norm. We use the se- 
quence of q + 1-tuples of sets 7(n) = (r°, . . . .V^^) to specify boundary conditions 
by imagining that all particles in F^ are of type j for j = 1, g and none occupies 
F°. This defines a restricted ensemble and the corresponding probability measure 
is denoted by /U„ = /Xn^"^''^'^. The choice of b.c.s 7(n) is understood to be fixed. 
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Let's consider the Ising-Potts model in an n x n lattice box B{n) with boundary 
conditions j € {1, 2, g} at a fixed inverse temperature /? with the corresponding 
probability measure Hg^^y It is well known that, as n — > oo, a unique, translation 

invariant infinite volume measure ij.^^'^ emerges as the weak limit of the sequence 

('^_B(n^)">i" '^^'^ define the order parameter = 6{p,q,d) as the excess density 
of the dominant color, the excess is measured from the symmetric value l/q. The 
model exhibits phase transition in the sense that in dimensions d >2 there exists a 
critical value < P^d) < oo such that for /3 < /?c, 0{p) := n^^''^ [ao = j] -l/q = 
but for /? > (3c, 0{(3) > 0, i.e., when the interaction becomes strong enough, the 
influence of the (arbitrarily) far away boundary still propagates all the way through 
the inside of the volume and creates a majority of j'-type particles. (Note that in 
the Ising model (q = 2), the spontaneous magnetization m* is equal to 9.) The 
probability measures fJ,^^'^ , j = l,2,...,q, describe in mathematical terms what 
we call 'pure' phases and which correspond to the saturated solutions in the initial 
example. 

Having chosen our model, let us formulate the goals of the analysis. Clearly, 
the main goal is to verify the predictions of the phenomenological theory, namely, 
that on the macroscopic scale the phases will be arranged according to some solution 
of the variational principle corresponding to minimal surface energy. First, however, 
the participating phases have to be found and identified. Moreover, as pointed out 
earlier, the previous statement contains a couple of implicit assumptions, such as the 
existence of the surface tension, the absence of transitional states (where one phase 
would smoothly go over into another one) the regularity of the interface boundaries, 
etc., all of which have to be justified from a microscopic point of view. 

1.3. Connection to minimal surfaces 

In this section a partially informal discussion of some examples will be pre- 
sented with the aim to make the close relation to minimal surfaces transparent. 





figure 1 

Consider the Potts model with q = 6 colors (states) in a three dimensional box 
with boundary condition i on the i-th face of the box. Naively, one might expect 
that all phases will try to occupy the region closest to the corresponding piece of 
the boundary, which would lead to a phase partition consisting of symmetric and 
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pyramid-like regions, as can be seen in figure 1, left. However, at least in the case 
when the surface tension is isotropic (which is presumably the case in the limit 
T t Tc), there exists a better configuration with lower total surface free energy. 
Recall that in this case our desired interface is simply a minimal surface spanned by 
the edges of the box. A picture of the well known solution to this problem can be 
seen in figure 1, right. In order to be able to discuss this example at temperatures 
< T < Tc, we have to make certain assumptions about the surface tension t. We 
assume that the sharp simplex inequality holds, that the value of r is minimal in 
axis directions and that r increases as the normal vector moves from say (0, 0, 1) 
to (1, 1, 1). (Although these assumptions are very plausible, none of them has been 
proved in dimensions d>2>). Under these hypotheses, we conjecture that the phase 
partition at moderate subcritical temperatures looks like in figure 2, left. In the 
limit T I 0, only two phases survive, as shown in figure 2, right. At T = 0, there 
is no reason for the middle plane to stay centered, in fact, any horizontal plane is 
equally likely. 





figure 2 

In the next example we consider the three dimensional Ising model with free 
boundary conditions below T^, conditioned on the event that the average magneti- 
zation is positive and does not exceed m* — e, where £ is a sufficiently small positive 
number and m* denotes the spontaneous magnetization. It is natural to conjecture 
that a minimizer of the corresponding variational problem is a droplet attached 
symmetrically to one of the corners of the box. 

The single bubble sitting in one of the corners is filled with the minus phase 
and in the rest of tlic^ box we sec the plus phase. The size of the bubble is determined 
by e and its internal boundary coincides with the corresponding piece of the surface 
of the Wulff crystal. 

Another Wulff-type problem arises by conditioning the (/-states Potts model 
(with say 9 > 4) to have a moderate excess of colors 2 and 3 while imposing 1- 
boundary conditions on the entire box. In this case it is conceivable that a so-called 
"double bubble" is created, consisting of two adjacent macroscopic droplets filled 
with the (pure) phases 2 and 3, respectively. The double bubble is swimming in 
the phase 1 which fills the rest of the box. Of course, we might have an excess of 
color 4 as well; in this case a further bubble will presumably appear which will be 
attached to the previous two bubbles. 

For related variational questions concerning soap films and immiscible fluids, 
see [30]. 
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In fact, by studying questions concerning phase boundaries we are very quickly 
confronted with the theory of minimal surfaces, such as the Plateau problem, cor- 
responding to anisotropic surface measures. Let O be a bounded open set in 
with smooth boundary and lot 7 bo a Jordan curve drawn on dVl which separates 
d fl into two disjoint relatively open sets r+ and T~ . Typical configurations in the 
Ising model on a fine grid in fl with plus b.c.s on r+ and minus b.c.s on F" will 
exhibit two phases separated with an interface close to a minimal surface which is 
a global solution of the following Plateau type problem: 



minimize / t{us{x)) dlif' ^{x) : 5 is a surface in O spanned by 7 
Js 

where i/s(x) is the normal vector to S at x. We remark that it is conjectured that, 
as the temperature approaches Tc from below, the surface tension r becomes more 
and more isotropic and it is conceivable that the solution of the above minimization 
problem approaches the solution of the classical (isotropic) Plateau problem. 

1.4. Further background 

There is a beautiful and extremely useful way to decompose the Ising-Potts 
model into a certain bond percolation model, called FK-percolation and some simple 
'coloring' procedure discovered by Fortuin and Kasteleyn [21]. Consider the g-state 
Potts model with mixed boundary conditions 7 in a finite lattice box B at inverse 
temperature /3 and set p = 1 — e~^. Consider a bond percolation model (called FK 
percolation) specified by the following formula: 



where 77 is a bond configuration with the property that there is no open connec- 
tion between differently colored boundary parts, Pp is the usual Bernoulli measure 
with parameter p, #(77) denotes the number of clusters (conn, components) in the 
configuration rj with the rule that identically colored boundary parts (and their con- 
nected components) count as one single cluster. Finally Zg'^'"^ is the appropriate 
normalizing constant. 

In the second step we assign colors to every cluster (and their sites) as follows: 
the boundary pieces inherit the color of the boundary condition, the remaining 
clusters will get one of the colors 1,2, ...,q with probability 1/q each independently 
from each other. The distribution of the coloring of the sites corresponds exactly 
to the Potts model. Note that in the case of free or constant b.c.s. there is no 
constraint prohibiting open connections, and indeed these measures, denoted by 
^b'^ '' (free b.c.s.) and <j)™'''''? ("wired" b.c.s.) behave very similar to regular 
Bernoulli percolation. Their thermodynamic limits and exist as the 

box size tends to infinity and we can define the percolation probability as usual 
^*(p) = ^ 00], for * free or wired. It is easy to check that the order 

parameter 6{P) of the Ising-Potts model agrees with 6'^" {(3), and correspondingly 
the FK model exhibits a percolation phase transition. Further it is known [23] , that 
6'^ {(3) = 0^ {(}) for all but at most countably many values of (i and it is conjectured 
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that the equality is vahd for all values except possibly the critical point 0c = Pcil)- 
Moreover this condition is equivalent to the equality of the thermodynamical limits: 
^ w,/3,9 _ ^Ll3,q_ ■\Yg define the set of "regular" inverse temperatures by 

ll{q) = {(3>0;e^-{P)^9f{P)}. 

Although the status of the bonds are dependent, their correlation tends fast to 
zero as the distance between them becomes large and this holds for both the sub 
and supercritical phase of FK percolation. This property is of crucial importance 
in our large deviation analysis since in the original Potts model there exists no 
corresponding asymptotic independence when /3 > /3c. 

Our results arc valid above the so called slab-threshold (3c = (3c{q, d), introduced 
in [34]. This threshold is conjectured to agree with the critical point and at least in 
the case of percolation (q= 1) this have been proved by Grimmett and Marstrand 
[24]. It is possible to characterize this threshold as the smallest value such that 
when /3 exceeds it, it is possible to find a > and L > 1 such that at least in the 
center of the of slabs S{L,n) = [—L,L] x [—n,n]'^~^ fl Z** there is "uniform long 
range order", i.e., 

inf inf x y > 0. 

n>l x,yeS{L,an) ^(^-'"■> 1 "J 



It has been proved in [34] that above (3c, a always can be chosen to be one, guar- 
anteeing a strictly positive probability (uniformly in n) for connections within a 
sufficiently (depending on (3) thick slab. This property is crucial for establishing 
the basic properties of supercritical FK percolation; the existence of a unique cross- 
ing cluster in a box, its omnipresence, the concentration of its density around the 
percolation probability, the exponential tail decay of the diameter of other clus- 
ters in the box, etc. These properties are then used to establish a renormalization 
scheme which is essential for the large deviation analysis of this and the Ising- Potts 
model. 



2. The results 



2.1. Historical remcirks 

Before we start with the presentation of our results we give a brief summary of 
the previous work on this subject. As we have already mentioned, large deviations 
theory plays an important role in this context and not surprisingly the first efforts 
were devoted to the study of large deviations of the empirical magnetization in the 
Ising model, i.e., the average value of the spins in a large box. A volum,e order 
large deviation principle (LDP) has been established for the Ising model by various 
authors: Comets, Ellis, Follmer, Orey, OUa [12, 15, 17, 32]. The corresponding rate 
function has been found to vanish in [— to*,to*] where to* denotes the spontaneous 
magnetization. In fact, it was suspected that the correct order of decay is expo- 
nential to surface order. Indeed, Schonmann [35] found a proof of this conjecture, 
valid for any dimensions and low enough temperatures and Chayes, Chayes and 
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Schonmann extended the result for the supercritical /3 > [3,, regime in the two di- 
mensional case. Follmer and Ort [18] investigated this phenomenon on the level of 
empirical measures. Finally, inspired by the work of Kesten and Zhang [29] on re- 
lated questions in percolation, Pisztora [34] established surface order upper bounds 
for the remaining dimensions d>2> above the slab-threshold j3c, introduced in the 
same work, which is conjectured to agree with the critical point (3c- In that work 
a renormalization scheme has been developed for supercritical Fortuin-Kasteleyn 
percolation (or random cluster model) in conjunction with a stochastic domination 
argument (generalized and improved in [31]) which allows to control the renormal- 
ized process, and so, the original one. 

The monograph of Dobrushin, Kotecky and Shlosman [14] opened the way to 
the rigorous study of the phase separation phenomenon creating thereby an immense 
interest and activity which lasts up to the present time. Their analysis, which 
provided the first mathematical proof of phase separation, had been performed in 
the context of the Ising model. The main tool of their work is the cluster expansion, 
which, on the one hand allowed the derivation of results much finer than necessary 
to verify the Wulff construction, on the other hand it restricted the validity of the 
results to two dimensions and low temperatures. Significant improvements of these 
results in two dimensions have been derived by Pfister [33] , Alexander, Chayes and 
Chayes [4] (treating percolation), Alexander [3], loffe [26, 27]. Finally loffe and 
Schonmann [28] extended the results of [14] up to Tc- 

The next challenge was to analyze phase separation for short range models in 
higher dimensions. The additional difficulties came mainly from two sources. First, 
new techniques have to be developed to avoid the use of perturbative methods (such 
as the cluster expansion) which severely limit the applicability of the arguments and 
methods which are specific to two dimensions only (duality). Second, the emerging 
geometry is far more complex than in two dimensions and this requires the use of 
new tools and ideas. The complexity of the geometry causes problems also within 
the probabilistic analysis (for instance the lack of the skeleton technique for surfaces) 
and even the correct formulation of the results is far from obvious ("hairs"). 

The first issue has been resolved by the application of the aforementioned 
renormalization technology from [34]. Renormalization arguments lie at the heart 
of the proof of much of the intermediate steps (for instance exponential tightness, 
decoupling) and even in the remaining parts they play an important role usually in 
combination with geometric arguments (interface lemma, etc.). 

To handle the geometric difficulties, the use of appropriate tools from geometric 
measure theory has been introduced in the works of Alberti, Bellettini, Bodineau, 
Butta, Cassandro, Presutti [2, 5, 6] and by Cerf [9]. In these works also a novel and 
very general large deviation framework have been proposed to tackle the problem. In 
fact, this framework turned out to be crucial for the success of the entire approach. 
It is the work of Cerf [9] in which the first complete analysis of phase separation in 
a three dimensional model have been achieved, namely the asymptotic analysis of 
the shape of a large finite cluster (Wulff problem) in percolation. 

The results presented in this article have been derived in the works [34, 9, 10, 
11]. It should be mentioned that in an independent work [7] Bodineau carried out 
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an analysis of the Wulff problem in the Ising model with conclusions slightly weaker 
than the results appearing in [10]. 

Finally, for current developments in the field we refer the reader to the preprint 
[8] and the references therein. 

2.2. Statement of the results 

Range of validity of the results. Our results for the Ising-Potts models hold in 
the region: d>3, q€N\{0,l}, f3 > Xiq, d), f3 € U{q, d). 

At this point it is natural to comment on the case of two dimensions. Although 
most of our results should hold for d = 2, there are several points in the proofs which 
would require a significant change, making the proofs even longer. The main reason, 
however, for not to treat the two dimensional case is that the natural topology for 
the LDP-s in d = 2 is not the one we use (which is based on the distance distLi) 
but a topology based on the Hausdorff distance. 

Surface tension. From FK percolation we can extract a direction dependent 
surface tension t{v) = T{p,q,d,i'), cf [10]. For a unit vector v, let A be a unit 
hypersquare orthogonal to v, let cyl A be the cylinder A + Riv, then t{v) is equal 
to the limit 



The function r satisfies the weak simplex inequality, is continuous, uniformly bounded 
away from zero and infinity and invariant under the isometrics which leave in- 
variant (see section 4 in [10] for details). 

Identification of the phases. The typical picture which emerges from the Potts 
model with mixed b.c.s. at the macroscopic level is a partition of O in maximal q 
phases corresponding to the dominant color in that phase. The individual phases 
need not be connected. In order to identify the phases we choose first a sequence 
of test events which we regard as characteristic for that phase. More specificly, 
for j = 1,2, ...,(7 we select events En' defined on a n x n lattice box such that 
/xi^'^''^'' [eI^^] ^ 1 as n ^ oo. We may also assume that E^^^ n E^J"* = for j i. 
For instance, one natural choice is to require that the densities of the different colors 
in the box do not deviate more than some small fraction from their expected value. 
This will guarantee that the right mixture of colors occurs which is typical for that 
particular pure phase. Alternately, we may request that the empirical measure 
defined by the given configuration is close to the restriction of jJi^'^''^ to A(n) with 
respect some appropriate distance between probability measures, etc. 

For a; e M** and r > we define the box A(a;, r) by 



lim 

n— »cx> TV 



1 



log$Si« 



/inside n cyl A there exists a finite set of closed edges E ' 
cutting ncylA in at least 2 unbounded components and 
the edges of E at distance less than 2d from the boundary 

Vof n cyl A are at distance less than 2d from nA , 



K{x,r) = { 




r/2 <yi- Xi < r/2, i = 1,. . . ,d 



and we introduce an intermediate length scale represented by a fixed function / 
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N — > N satisfying 

lim n/ f{nY~^ = lim /(n)/ logn = oo. (1) 

n — ^oo n — *oo 

Given a configuration in f2„, we say that the point a; G belongs to the phase j, if 
the event En^ occurs in the box A(a;, f{n)/n). For j = 1, 2, .., q, we denote by An^ 
the set of points in O belonging to the phase j and set An^ = O \ UjAn^ (indefinite 
phase) . 

The random partition of il, A„ = (A" , A^, . . . , ^^), is called the empirical 
phase partition. Our first result shows that up to super-surface order large devia- 
tions, the region of indefinite phase A^ has negligible density, i.e., the pure phases 
fill out the entire volume. 

Theorem 2.1 Letd>3,q£n\ {0, 1}, P > X, P & d). For 6 > Q, 
lim sup log /i„[vol(A°) > ^] = -00. 

n — >oo ^ 

Although Theorem 2.1 guarantees that the pure phases fill out the entire vol- 
ume (up to negligible density) but it does not exclude the possibility that the 
connected components of the pure phases arc very small. For instance, they could 
have a diameter not much larger than our fixed intermediate scale (in which case 
they would be invisible on a macroscopic scale.) If this happened, the total area 
of the phase boundaries would be exceedingly high. Before stating our next result, 
which will exclude this possibility, we introduce some geometric tools. 

We define a (pseudo) metric, denoted by distLi, on the set B(f2) of the Borel 
subsets of by setting 

VAi,^2eB(0) distLi(Ai,A2) = vol(AiAA2). (2) 

We consider then the space of phase partitions P{^, q) consisting of g -|- 1-tuples 
{A^ ,A^, . . . , A'i) of Borel subsets of f2 forming a partition of O. We endow P(fi, q) 
with the following metric: 

distp((^o,...,A«),(BO,...,B«)) = distLi(^%B^). 

i=0,...,9 

In order to define the surface energy X of a phase partition An, we recall some 
notions and facts from the theory of sets of finite perimeter, introduced initially 
by Caccioppoli and subsequently developed by De Giorgi, see for instance [22, 16]. 
The perimeter of a Borel set E of is defined as 

V{E) = sup { J^divf{x)dx : f e Co~(M^5(l))} 

where C^{R'^, B{1)) is the set of the compactly supported C°° vector functions 
from R'^ to the unit ball B{1) and div is the usual divergence operator. The set E is 
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of finite perimeter if V{E) is finite. A unit vector u is called the measure theoretic 
exterior normal to E at x if 



Let E he a set of finite perimeter. Then there exists a certain subset of the topo- 
logical boundary of E, called the reduced boundary, denoted by d*E, with the same 
d — 1 dimensional Hausdorff measure as dE, such that at each x G d*E there is a 
measure theoretic exterior normal to E at x. For practical (measure theoretic) pur- 
poses, the reduced boundary represents the boundary of any set of finite perimeter, 
for instance, the following generalization of Gauss Theorem holds: For any vector 



(For more on this see e.g. the appendix in [10] and the references there.) 

The surface energy I of a phase partition {A^A\, . . . , ^^) e P(0, q) is defined 
as follows: 

- for any (A°, A\ . . . , A^) such that either A° 7^ or one set among A\ . . . , ^« 
has not finite perimeter, we set X{A°, . . . , A') = 00, 

- for any (A°, A\ . . . , A«) with A° = and A^,...,A^ having finite perimeter, 



Note that X depends on r and the boundary conditions 7 = (F"'^, F^). The 
first term in the above formula corresponds to the interfaces present in Q, while 
the second term corresponds to the interfaces between the elements of the phase 
partition and the boundary F. It is natural to define the 'perimeter of the phase 
partition by using the same formula with t replaced by constant one. Since t is 
uniformly bounded away from zero and infinity at any temperature, the surface 
energy can be bounded by a multiple of the perimeter and vice versa. It is known 
that the surface energy I and the perimeter V are lower semi continuous and their 
level sets of the form T~^[0, K] are compact on the space (S(M'^), distLi). 

The next result states that up to surface order large deviations (and the con- 
stant can be made arbitrarily large by adjusting the bound K below) the empirical 
phase partition will be (arbitrarily) close to the set of phase partitions with perime- 
ter not exceeding K. 




function / in (»'',»'') 




we set 




Theorem 2.2 Letd>3, qeN\ {0, 1}, (3 > l3c, l3 e U{q, d). For ^ > 0, 



lim sup 



1 



log/z„ distp(A„, X-^[0,iV']) > (5 <-cK. 



n' 



d-l 
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Our fundamental result is a large deviation principle (LDP) for the empirical 

phase partition , Al^, . . . , A"^). 

Theorem 2.3 The sequence {^n)^^^ — ((^n: ■ • ■ ; ^n))„gN °f empirical 
phase partitions of Ct satisfies a LDP in {P{fl,q), distp) with respect to jin with 
speedn'^~^ and rate function I— vimip^Q^q) I, i.e., for any Borel subsetE of P{fl,q), 



■ inf T + min X < lim inf —t—t log u„ 

o pro n^ n—^oo n"--^ 



P{n,q) 

E 



< lim sup --^ log fin 



n—^oo ^ 



< — inf X + min X. 



Note that the constant minp(o^g)X is always finite. Every large deviation re- 
sult includes a (weak) law of large numbers; here the corresponding statement is 
as follows: Asymptotically, the empirical phase partition will be concentrated in 
an (arbitrarily) small neighborhood of the set of partitions minimizing the surface 
energy. In other words, on the macroscopic level, the typical phase partition will 
coincide with some of the minimizers of the variational problem, in agreement with 
the phenomenological prediction. Of course, the LDP states much more than this, 
in particular we will be able to extract similar statements for restricted ensembles. 
Recall that imposing mixed boundary conditions is not the only way to force the sys- 
tem to exhibit coexisting phases. In the Wulff problem in the Ising model context, 
for instance, a restricted ensemble is studied which is characterized by an artifi- 
cial excess of say minus spins in the plus phase. Technically this can be achieved 
by conditioning the system to have a magnetization larger than the spontaneous 
magnetization while imposing plus b.c.s. 

The next result describes the large deviation behavior of the phase partition 
in a large class of restricted ensembles. Although it is a rather straightforward gen- 
eralization of Theorem 2.3, we state it separately because of its physical relevance. 

Let {Gn)n>i be a sequence of events, i.e., sets of spin configurations, satisfying 
the following two conditions: first there exists a Borel subset G of P{^, q) such that 
the sequence of events {Gn)nen and {{An G G})„£n are exponentially equivalent, 
i.e., 



limsup— T-j- log Hn 



Gn A {An G G} = -oo, (3) 



where A denotes the symmetric difference. Second, the following limit exists and 
is finite: 

Xg = lim — ^ log fin[Gn] > -00. (4) 
n — *oo Ji"- ^ 

The sequence of events (Gn)n>i determines a restricted (conditional) ensemble. 
Note that if 

inf X = inf X > — oo, (5) 
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then Theorem 2.3 imphes that (4) is satisfied, with Ig = inf I. 

Theorem 2.4 Assume that the sequence {Gn)n>i satisfies (3) and (4) and define 
for each n> 1 the conditional measures 

Then the sequence of the empirical phase partitions offl satisfies a LDP 

in {P{fl,q), distp) with respect to fi^ with speed n'^~^ and rate function I — Iq, 
i.e., for any Borel subset E of P{^l,q), 



- inf I + Jg < lim inf — 7— r log 

EnG " 

1 e 

< hmsup-^ log 11^ 

n— i-oo 

< - inf J + Jg- 

EnG 



i*™ ge 



Theorem 2.4 gives a rigorous verification of the basic assumption underlying 

the phenomenological theory, namely, that in a given ensemble, the typical con- 
figurations are those minimizing the surface free energy. A general compactness 
argument implies the existence of at least one such minimizer. However, in most 
examples one cannot say much about the minimizers themselves. (One notable ex- 
ception is the Wulff problem.) The difficulty stems from the fact that the surface 
tension r is anisotropic and almost no quantitative information about its magnitude 
is available. Moreover, the corresponding variational problems are extremely hard 
even in the isotropic case and the (few) resolved questions represent the state of 
the art in the calculus of variations. For instance, a famous conjecture related to 
the symmetric double-bubble in the three dimensional case with isotropic surface 
energy (perimeter) has only been resolved recently [25] and the asymmetric case 
remains unresolved (even in the isotropic case). 

We show next how Theorem 2.4 can be applied to the Wulff and multiple 
bubble problem. We take pure boundary conditions with color 1, that is, = F, 
F^ = • • • = F'' =0. Let S2, • • • , be g — 1 real numbers larger than or equal to 
{l-e)/q. We set 

Vi e {2, . . . , g} Vi= vol {Q.)e-^{si - (1 - e)/q). 

We define next the events 

VneN Gn = {ii&{2,...,q} 5„(i) > 

and the collection of phase partitions 

<G{v2, ...,Vq) = {A={Aa,Ai,...,Aq)e P{n, q) : vol (A2) > t-a, . . . , vol (A,) > }. 
It can be shown that the sequences of events 

(G„)„eN and {A^ G G(u2, . . . , Vq))nen 
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arc exponentially equivalent, i.e., they satisfy the condition (3). In order to en- 
sure condition (5), we suppose that the minimum of the surface energy X over 
G(v2, . . . , Vq) is reached with a phase partition having no interfaces on the bound- 
ary r. More precisely, wc suppose that the following assumption is fulfilled. 
Assumption. The region and the real numbers vit- - ^Vq are such that there 
exists A* = {Aq, Al,..., A*) in G(w2, ■ ■ ■,Vq) such that 

1(1*) = mm{I{A);AeG{v2,...,Vg)} 
Wie{2,...,q} d2iA*,r)>o. 

We expect that this assumption is fulfilled provided the real numbers V2,---,Vq 
are sufficiently small (or equivalently, S2, ■ ■ ■ ,Sq are sufficiently close to (1 — 0)/q), 
depending on the region This is for instance the case when q = 2. Indeed, let 
Wt be the Wulff crystal associated to r. We know that Wr is, up to dilatations and 
translations, the unique solution to the anisotropic isoperimetric problem associated 
to T. For V2 sufficiently small, a dilated Wulff crystal xo+aoWr of volume V2 fits into 
f2 without touching F , and the phase partition A* = (0 , \ (xq + ccq Wt ),xo + ao Wt ) 
answers the problem. In the case q > 2,we expect that a minimizing phase partition 
corresponds to a multiple bubble having q — I components. 

Under the above assumption, we claim that the collection of phase partitions 
G{v2, . . . ,Vq) satisfies (5). For A > 1, we define 

A*{X) = (9,n\ U XA*,XA*2,...,XA;). 

2<i<q 

Since by hypothesis the sets A^,. . . ,A* are at positive distance from F, for A larger 
than 1 and sufficiently close to 1, the phase partition A*{X) satisfies 

A*{X) e G{X''v2,...,X\) C G(v2,---,vq) 

and moreover X{A*{Xj) = X'^~^X{A*). Sending A to 1, and remarking that 
G{v2, ■ ■ ■ ,Vq) is closed, we see that G{v2, ■ ■ ■ ,Vq) satisfies (5). Thus we can ap- 
ply Theorem 2.4 with the sequence of events {Gn)neN, thereby obtaining a LDP 
and a weak law of large numbers for the conditional measures fj.^ = /Lt„(- |G„). 
In the particular case q = 2, we obtain again the main result (Wulff problem) of 
our previous paper [10]. In the more challenging situations q > 2, the imresolvcd 
questions concerning the macroscopic behavior of such systems belong to the realm 
of the calculus of variations. 
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